9ed9b1e6-42cf-4e8c-ac21-4b621403ef9b 59 6 5 auto

&[DATE] &[TIME] - &[FILENAME]

CWD C:\SmathFile CurrentDirectory :

M 1 -^br * ecal M 1 -^br * 11 line ≡

The [lpc] Spline Project

Part 1. get the 's' spline coefficients

t0 1 time :

L H N xd range discretise U 0 : dx H L - N / : i 1 N range U i el L dx - (dx i * + : for U 21 line 41 line :

L 0 :

H 10 :

Num 100 :

U L H Num xd :

n U rows :

n 100

x φ 11 x 5 - (2^+ / :

Data i 1 n range vx i el U i el : vy i el U i el φ : 21 line for vx vy augment 21 line :

Data

X Data 1 col :

Y Data 2 col :

N n 1 - :

N 99

r 2 N range :

Prep.0 Install the constant vector 'br'

i r br i el Y i 1 + el Y i el - X i 1 + el X i el - / (Y i el Y i 1 - el - X i el X i 1 - el - / (- : br N 1 + el 0 : 21 line for

Prep.1 Install the tridiagonal matrix

M 11 el 1 :

M N 1 + N 1 + el 1 :

0. Construct the tridiagonal linear spline matrix [lspline]

i r M i i 1 - el X i el X i 1 - el - 6 / : for

i r M i i el X i 1 + el X i 1 - el - 3 / : for

i r M i i 1 + el X i 1 + el X i el - 6 / : for

1. Construct the tridiagonal parabolic spline matrix [pspline]

M 12 el 1 - :

M N 1 + N el 1 - :

2. Construct the tridiagonal cubic spline matrix [cspline]

M 12 el X 1 el X 3 el - (X 3 el X 2 el - / (:

M N 1 + N el X N el X N 2 - el - (X N 2 - el X N 1 - el - / :

M 13 el X 2 el X 1 el - (X 3 el X 2 el - (/ :

M N 1 + N 1 - el X N 1 - el X N el - (X N 1 - el X N el - (/ :

M 1 -^br * eval solve the tridiagonal matrix ≡

Spool to file the spline coefficcients M 1 -^br * eval CWD spline wfile 1

Here, we collect the 'cspline'. Disable the 'cspline' tridiagonal to collect 'pspline' ... disable 'pspline' tridiagonal to collect 'lspline.

s CWD spline rfile :

s rows 100

Part 2. construct the cumulative integral

vx X :

vy Y :

s rows 100

Prepare the pseudo-integrand s vx vy fintegrate i 1 n 1 - range h i el vx i 1 + el vx i el - : f i el h i el 2 / vy i el vy i 1 + el + (* h i el (3^24 / s i el s i 1 + el + (* - : 21 line for f 21 line :

f s vx vy fintegrate :

v vx 1 f rows 11 submatrix :

f rows v rows 12 mat 999912 mat

view the pseudo-integrand

v f augment

Export the cumulative integral Σcumul n vx rows 1 - : Sum 1 el 0 : Sum n 1 + el 0 : i 1 n range Sum i 1 + el Sum i el f i el + : for vx Sum augment 51 line :

x MyInt x φ x int maple : 5 - x + atan

sl_p_c spline cumulative integral

Σcumul 0 2 1 sys

Symbolic cumulative Integral

x 0 ≥ ( x 10 ≤ ( & x MyInt 0 MyInt - Maple Symbolic if 0 2 1 sys

Part 3. The complimentary scalar dspline

s X Y t dspline t X 2 el ≤ k 1 : t X X rows 1 - el ≥ k X rows 1 - : k 1 : t X k 1 + el > k k 1 + : while 21 line if if c Y k 1 + el Y k el - X k 1 + el X k el - / s k 1 + el 2 s k el * + (X k 1 + el X k el - 6 / (* - : s k 1 + el s k el - X k 1 + el X k el - / (t X k el - (2^2 / * s k el t X k el - (* + c + 31 line :

x D s X Y x dspline :

x Dsymbolic x φ x 2 diff :

x D2 x D x diff :

<= Smath does NOT sub-derive

The scalar dspline

0 x ≤ ( x H ≤ ( & x D scalar dspline if 0 2 1 sys

2nd derivative of [lpc] spline

X s augment 0 2 1 sys

Symbolic derivative

0 x ≤ ( x H ≤ ( & ( x Dsymbolic D2 Symbolic if 0 2 1 sys

1 time t0 - 25.4 sec *

Project Note: In this demo, we have used a function. In real life you will just have a data set and no 2nd symbolic unless you have fitted by some method. In the case of "just a data set", both derivatives will be 0K. Neither 1rst/2nd derivatives need smoothing.

Part 4. The [lpc] spline interpolation

The cubic spline [lpc], scalar interpolation, plot scalar s X Y t InterpSpline t X 1 el ≤ idx 1221 mat : t X N el ≥ idx N N 1 + 21 mat : i 1 : X i el t < i i 1 + : while idx i 1 - i 21 mat : 31 line if if hx X idx 2 el el X idx 1 el el - (: hy Y idx 2 el el Y idx 1 el el - (: hs s idx 2 el el s idx 1 el el - (: 31 line a hs 6 hx * / : b s idx 1 el el 2 / : c hy hx / hx 6 / s idx 2 el el 2 s idx 1 el el * + (* - (: d Y idx 1 el el : 41 line a t X idx 1 el el - (3^* b t X idx 1 el el - (2^* + c t X idx 1 el el - (* + d + 41 line :

x Spline s X Y x InterpSpline :

"Proofing" the cubic spline project

0 x ≤ ( x H ≤ ( & ( x Spline Job done ! if

Sources and credits: 1. Mathsoft for the tridiagonal systems. 2. Robert Adair for the scalar algos: InterpSpline(s,X,Y,t) ... dspline(s,X,Y,t) ... fintegrate(s,vx,vy)

jmGiraud: 2016_02_18

x Deriv x Spline (x diff :

Conclusion:

Smath does not recognise Deriv(x) as a scalar function. The Smath D(x) operator sets some algebraic rules book style. The Ridder/Brent operator is missing.